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Dynamics of an Ion Coupled to a Parametric Superconducting Circuit 
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Superconducting circuits and trapped ions are promising architectures for quantum information 
processing. However, the natural frequencies for controlling these systems - radio frequency ion 
control and microwave domain superconducting qubit control - make direct Hamiltonian interactions 
between them weak. In this paper we describe a technique for coupling a trapped ion’s motion to 
the fundamental mode of a superconducting circuit, by applying to the circuit a carefully modulated 
external magnetic flux. In conjunction with a non-linear element (Josephson junction), this gives 
the circuit an effective time-dependent inductance. We then show how to tune the external flux to 
generate a resonant coupling between the circuit and ion’s motional mode, and discuss the limitations 
of this approach compared to using a time-dependent capacitance. 
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1. INTRODUCTION 

Superconducting circuits and trapped ions have dis¬ 
tinct advantages in quantum information processing. 
Circuits are known for fast gate times, flexible fabrica¬ 
tion methods, and macroscopic sizes, allowing multiple 
applications in quantum information science [1-4]. Un¬ 
fortunately, they have short coherence times and their 
decoherence mechanisms are hard to address [5]. Trapped 
ions, on the other hand, serve as ideal quantum mem¬ 
ories. Indeed, the hyperfine transition displays coher¬ 
ence times on the order of seconds to minutes [6-9], while 
high fidelity state readout is available through fluores¬ 
cence spectroscopy [10, 11]. Unfortunately, ions depend 
mainly on motional gates for interactions [12-15]; these 
are correspondingly slow, susceptible to motional heat¬ 
ing associated with traps[16-18], and occur only at a 
short - dipolar - range. The distinct advantages of ion 
and superconducting systems therefore motivates a hy¬ 
brid system comprising both architectures, producing a 
long-range coupling between high-quality quantum mem¬ 
ories. 

Early proposals for hybrid atomic and solid state sys¬ 
tems [19, 20] have yet to be implemented experimentally. 
Other approaches involve coupling solid state systems 
to atomic systems with large dipole moments, such as 
ensembles of polar molecules[21] or Rydberg atoms[22]. 
Although the motional dipole couplings with the elec¬ 
tric field of superconducting circuits can be several hun¬ 
dred kHz, these systems suffer from a large mismatch 
between motional (~ MHz) and circuit (~ GHz) fre¬ 
quencies. This causes the normal modes of the coupled 
systems to be either predominantly motional or photonic 
in nature, thereby limiting the rate at which information 
is carried between them. Implementation of a practical 
hybrid device therefore requires something additional. 

Parametric processes allow for efficient conversion of 
excitations between off-resonant systems. In the field of 
quantum optics, they are widely used in the frequency 
conversion of photons using nonlinear media [23, 24]. In 
the realm of superconducting quantum devices, paramet¬ 


ric amplifiers provide highly sensitive, continuous readout 
measurements while adding little noise [25-28]. Para¬ 
metric processes can also been used to generate control¬ 
lable interactions between superconducting qubits and 
microwave resonators [29-31]. In the context of hybrid 
systems, Ref. [32] presents a parametric coupling scheme 
between the resonant modes of an LC circuit and trapped 
ion. The ion, confined in a trap with frequency uJi, is 
coupled to the driven sidebands of a high quality factor 
parametric LC circuit whose capacitance is modulated at 
frequency Wiq —This gives rise to a coupling strength 
on the order of tens of kHz for typical ion chip-trap pa¬ 
rameters. A different approach proposed in Ref. [33] is 
based on a position-charge interaction that is quadratic 
in the charged particle’s position, so that driving of its 
motion produces a parametric coupling. In this approach 
the circuit couples to an electron rather than an ion, lead¬ 
ing to an enhanced coupling strength (~ 1 MHz) due to 
its reduced mass. 

In this manuscript we describe an alternative paramet¬ 
ric driving scheme to produce coherent interactions be¬ 
tween atoms and circuits. By using a time-dependent ex¬ 
ternal flux, we drive a superconducting loop containing 
a Josephson junction, which causes the superconductor 
to act as a parametric oscillator with a tunable induc¬ 
tance. By studying the characteristic, time-dependent 
excitations of this system, we show how to produce a res¬ 
onant interaction between it and the motional mode of a 
capacitively coupled trapped ion. Although in principle 
our approach could be used to produce a strong coupling, 
we find that, in contrast with capacitive driving schemes, 
the mismatch between inductive driving and capacitive 
(charge-mediated) interaction causes a significant loss in 
coupling strength. 

The manuscript proceeds as follows: In the first sec¬ 
tion, we describe the experimental setup of the circuit 
and ion systems and motivate the method used to couple 
them. To understand how the Josephson non-linearity 
and external flux affect the circuit, we transform to a 
reference frame corresponding to the classical solution 
of its non-linear Hamiltonian. We then linearize about 
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FIG. 1: Left: An rf-Squid with Josephson energy Ej and 
junction capacitance Cj driven by a time-dependent external 
flux The outer loop with inductance L contributes an 

energy i?_L. Right: Schematic of an ion confined by a capacitor 
with capacitance C. The capacitor is connected in parallel to 
the rf-SQUID system. 


this solution, resulting in a Hamiltonian that describes 
fluctuations about the classical equations of motion and 
corresponds to an LC circuit with a sinusoidally varying 
inductance. This periodic, linear Hamiltonian is charac¬ 
terized by the quasi-periodic solutions to Mathieu’s equa¬ 
tion. From these functions we define the ‘quasi-energy’ 
annihilation operator for the system and transform to 
a second reference frame where it is time-independent. 
Finally, we derive the circuit-ion interaction in the in¬ 
teraction picture, which allows us to directly compute 
the effective coherent coupling strength between the sys¬ 
tems. We conclude by comparing our results to previous 
work[32] (a capacitive driving scheme) and analyzing why 
inductive modulation is generically ineffective for capac¬ 
itive couplings. 


2. PHYSICAL SYSTEM AND HAMILTONIAN 

We begin with a basic physical description of the ion- 
circuit system, and as in Ref. [34], construct the classical 
Lagrangian before deriving the quantized Hamiltonian. 
We consider a single ion placed close to two capacitive 
plates of a superconducting circuit. We assume that the 
ion is trapped in an effective harmonic potential with a 
characteristic frequency Wj. in the direction parallel to the 
electric field between the two plates[10]. The associated 
ion Lagrangian is 


Ao„(z,i,f) =, (1) 

where m is the ion mass and z its displacement from equi¬ 
librium. Since the LC circuit will only couple to the ion 
motion in the z direction, we ignore the Lagrangian terms 
associated with motion in the other axial directions. 


The circuit interacting with the ion contains a Joseph¬ 
son junction [35] (Fig. 1) shunted by an inductive outer 
loop with inductance L. This configuration is known as 
a radio frequency superconducting quantum interference 
device (rf-SQUID)[2]. The SQUID is connected in par¬ 
allel to a capacitor C, whose plates are coupled to the 
nearby trapped ion. An external, time-dependent mag¬ 
netic flux through the circuit is used to tune the 

characteristic frequencies of this system. In our proposed 
configuration, the SQUID Lagrangian can be written in 
terms of the node flux (jj as [34] 

+ Ejcos{(j)/^o) - , 

(2) 

where 4>o = h/(2e) is the reduced flux quantum, and 
Cy = C + Cj is the effective capacitance of the circuit 
including the capacitance Cj of the Josephson junction. 
The Josephson energy is defined in terms of the critical 
current of the junction, Ej = Ic4>o- The Josephson term 
adds a non-linearity to the system oc Ej which will, in 
conjunction with the external driving, allow for coherent 
transfer of excitations between the otherwise off-resonant 
LC circuit (typical frequency ~ 1 — 10 GHz) and ion mo¬ 
tional mode (~ 1 — 10 MHz). Finally, we include the 
interaction between the two systems, associated with the 
ion motion through the electric field between the two 
plates in the dipole approximation: 

Ci = -e(J^^^z, (3) 


where d is the distance between the two plates and ^ ^ 
0.25 a dimensionless factor associated with the capacitor 
geometry [32]. 

With the full Lagrangian C = Cq+ Cion + Cj we follow 
the standard prescription to define the canonical vari¬ 
ables conjugate to z and cj): 


dC 

Pz = =mz 

az 

( 4 ) 

I = Try = Ue0-4-z. 

( 5 ) 

d(^ d 


These are identified as the momentum of the ion in the z 
direction and the effective Cooper-pair charge difference 
across the Josephson junction, respectively. We use the 
canonical Legendre transformation E[ = qcj) + pzZ — C to 
define the Hamiltonian and quantize the system, giving 


H{t)=H,on+Hq{t)EHl, (6) 


where 


Hj(t) 

Hi 


Pz ,1 2 ~2 


2Cy 


- Ej COs((^/^o) 


dC, 


-qz, 


( 7 ) 

^(0-<^.(^))^(8) 

(9) 
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As the canonical charge q of equation (5) has a term 
proportional to z, the ion Hamiltonian Hion gains an 
extra potential term proportional to Accounting for 
this term, we define the dressed trap frequency, 


2 2 
w,- = w. 


(PCY,m ' 


( 10 ) 


The latter correction is small compared to it is on the 
order of 50kHz for the parameters suggested in Ref. [32] 
{d ~ 25^m, Ce « 50fF and m « 1.5 x 10“^®kg the 
atomic mass of Beryllium). Finally, we note the resulting 
commutation relations, {4>,q\ = [z,Pz\ = ih. 

To motivate the method for coherently coupling the 
ion and circuit, we consider the Hamiltonian in the in¬ 
teraction picture, i.e., in the frame rotating with respect 
to Hq{t) + Hion- In this frame the Hamiltonian takes the 
form. 


= + (II) 

where b = annihilation opera¬ 

tor associated with excitations of the ion motional mode, 
and q{t) is propagated according to the time-dependent 
Hamiltonian Hq{t). Observe that if we set both 4)x{t) 
and Ej to 0 in the definition of Hq{t), the charge q{t) 
will oscillate as a harmonic oscillator, 

m ^ - ate-o‘), ( 12 ) 

with d the circuit annihilation operator analogous to b, 
and ojQ the undriven frequency LC, 

u;o = 1/x/L^. (13) 

In order to coherently exchange excitations between the 
systems, Hint{t) should only contain beam splitter-like 
terms, oE and while suppressing the excitation non¬ 
conserving terms, a)E and db. Yet since ojq ^ ooi, the 
o&I terms oscillate at a frequency comparable to the b^ 
terms, both of which have negligible effect in the rotating 
wave approximation (RWA). Our goal is therefore the 
following: design Ej and (pxit) such that dE contains a 
time-independent component in the interaction picture, 
while db contains only oscillating terms that drop out in 
the RWA. 


3. LINEARIZATION OF THE PARAMETRIC 
OSCILLATOR 

Before we can determine the parameters (px (t) and Ej 
allowing for coherent exchange of excitations, we must 
first bring Hq(t) into a more agreeable form. To begin, 
we linearize the circuit Hamiltonian about the classical 


solution to its equations of motion. This makes it possi¬ 
ble (in the next section) to identify effective annihilation 
and creation operators for the LC system, and to analyze 
their spectrum. 

We linearize by displacing the charge and flux variables 
by time-dependent scalars, qc{t) and (pcit), through the 
unitary transformations 

{7i(t) = 

(14) 

Specifically, we consider the LC circuit in terms of the 
displaced states, 

^{t))=Ul{t)U\{t)\^{t)) , (15) 


whose equation of motion satisfies 




Hq = ululHqUxU2 - ihUlUl'^U2 
■fcCt ^^2 




(16) 


Given [^,< 7 ] = ih, we have that U\qiji = q -\- qdt) 
and U\(j)U 2 = p + Pcit), while clearly UlpUi = (j> and 
lj\qlJ 2 = q- The displaced state Hamiltonian is therefore 

Hq{t) = ^l±^^+Vq{P + <Po(t)) 

+ {(j) +(t)c{t))dtqc{t) - qdtpcit), (17) 
where we have defined the nonlinear potential, 

Vqii’) = -Ej cos{p/po) + , (18) 


with (po = h/{2e). We now Taylor expand Vq about (pc{t), 
and collect all first and second order terms in q and p, 
giving 


Hq{t) = q {qc{t)/Cs - dtpc) + P (Vqipc) + dtqc) 
^2 


-k 


q 


V''iPc):;, 


P^ + RiP)-. 


2 Ce 2 ! 
where we have dropped all scalar terms, and 

R{P) = Y.~Plv!;'^\Po) 

fc>3 


(19) 


( 20 ) 


represents all higher order terms in the Taylor expansion. 

To complete the linearization, the displacements qc and 
pc are chosen so that the first order terms in q and p 
vanish, and therefore Hq is quadratic to leading order: 


dtpc = qc/C^ 
dtqc = -V^iPc). 


( 21 ) 
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These relations are the solutions to the classical, driven 
2 

Hamiltonian, Hq = and can be substi¬ 

tuted into each other to give 

dt(l)c + ujQ{(j)c + /3^o sin((/)c/^o)) = , (22) 

where 

/3 = LEj/~^l = ^ (23) 

90 

represents the strength of the non-linearity and ojq = 
\I\JLCy. is the bare LC resonance frequency. Sub¬ 
stituting from ( 21 ) and computing V”{(l)c) = j(l + 
P cos{(j)c/4’o))y we can now express the circuit Hamilto¬ 
nian as 

^ ^ ^(1+/3cos((/)c(t)/<^o))^^ (24) 

Note that we have dropped the higher order terms R{^) 
of equation (20). Understanding when this is valid will 
require us to express the flux operator (j) in the interaction 
picture, which we carry out below. A complete analysis 
is deferred to the appendix. 

With equation (24) we have converted to the Hamilto¬ 
nian of a harmonic oscillator with time-dependent induc¬ 
tance. Although it is explicitly dependent on the classical 
solution we note that engineering this function is 

rather straightforward. Given a desired the exter¬ 

nal driving 4>x{t) needed to produce it is given explicitly 
by equation (22) (see Fig. 2). For simplicity we assume 
that ipdt) satisfies the relation 

p cos{(j)c{t)/po) = V cosicodt ), (25) 

where we assume 77 < /3 so that pc is well defined at 
all t. This corresponds to an LC circuit with inverse 
inductance L~^ modulated at amplitude r] and frequency 
OJd, 

^ ++ ■ (26) 

In the following section we will see how the design param¬ 
eters rj and iOd determine the evolution of this system. 

4. TIME-DEPENDENT QUANTUM 
HARMONIC OSCILLATOR 

Although we now have a quadratic Hamiltonian for 
the SQUID, since Hq{t) = + 2^(1 + r] cos{oJdt))(j>‘^ is 

time-dependent, it is not immediately clear how to define 
its associated annihilation operator. To do so, we follow 
the approach of Ref. [36] to obtain an effective opera¬ 
tor a retaining useful properties of the time-independent 
case. Specifically, it satisfies [a, = 1 and, in the ap¬ 

propriate reference frame, is explicitly time-independent. 



FIG. 2: External flux drive px{t)/pQ required for a sinusoidal 
modulation of circuit inductance. The flux px{t) is defined 
explicitly by equation ( 22 ), where the classical solution pc{t) 
satisfies P cos{pcit)/ po) = 77 cos(a;df). The plotted function 
corresponds to drive parameters uid = 0 . 999 cjo , rj = 0.06 = 
(3/4) X P, where P = Ej{p^/L)~^ represents the strength 
of the non-linearity compared to the bare linear inductance. 
The relative scaling r]/P <l\s chosen so that pc is a smooth 
function of time; the 77 —>■ P~ limit corresponds to a saw-tooth 
wave with undefined second derivative. 

Further, in this frame the Hamiltonian can be written 
as h{dt0{t)){(d a + 1 / 2 ), where 9{t) is the effective phase 
accumulated by a in the Heisenberg picture. This new 
Hamiltonian commutes with itself at all times, allow¬ 
ing us to directly transform to the interaction picture 
in the next section. We will then describe how to control 
the spectral properties of q{t) in the interaction Hamilto¬ 
nian ( 11 ), producing the desired time-independent beam 
splitter-like interaction between circuit and ion motion. 

Following Ref. [36], we use unitary transformations to 
change Hq{t) so that it is of the form ~ g{t){aq'^ + Sp'^), 
where a and S are explicitly time-independent. To do 
this, we first analyze the classical equation of motion for 
the flux p associated with Hqit), as derived from equation 
(26), 

dtf{t)+ujl{l + 'qcos{uJdt))f{t) =0. (27) 

Since the drive has period r = 27r/a;d, from Floquet’s 
theorem [37, 38] we can find a quasi-periodic solution / 
satisfying 

/(O) = 1, 

/(i + r) = (28) 

In order for the solution to be stable, the characteris¬ 
tic exponent /i must be real valued. This imposes con¬ 
straints on the parameters Wd and 77 , which we discuss 
later. As we shall see, the properties of the function / 
will be closely related to the spectrum of q{t) in the in¬ 
teraction picture. 

It is useful to express / in polar form, 

f{t) = r(t)e*^(‘) , (29) 

where r > 0 and 9 is real valued. Substituting this re¬ 
lation into the characteristic equation (27) determines 
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equations of motion for r and 9, 


The final transformed Hamiltonian is therefore 


d^r = {{dtOY - Wo(l + 77COs(a;dt))) r , 

o = rd'iet + 2{dtr){dte). 


We define the associated Wronskian, which sets a char¬ 
acteristic frequency scale for the evolution, 


Hlc 


2 C^r^ 


2 




dtO 

W 





(38) 


= Im (re“*®5t (re*^)) 

= r^dtO. (31) 

We note that ^dtW is equal to the second line of (30), 
so dtW = 0 and W is time-independent. 

The above definitions are key to the transformations 
bringing Hq(t) into the desired form. Our first transfor¬ 
mation is 


where the last line follows from equation (31). 

Equation (38) allows us to define the effective annihi¬ 
lation operator in the standard way, 


■ /Ce»'. , . I i . 

so that the Hamiltonian is also equal to 

HLc{t)=hidt9)ia^ a+ 1/2). 


(39) 


(40) 


Us = exp{ixit)$'^/h ), (32) 

where 

X(t)^^^ (33) 

is the real part of the effective admittance, ^ ^. Using 

[0, q] = ih, we note 

UlqUs = q + 2 x^ 

uUUs = 0 

-ihUldtUs = ('9tx)0^ (34) 

After some algebra, the Hamiltonian of Eq. (26) is trans¬ 
formed to tjlHq{t)Us — ihUldtUs and becomes 

+ -^ (^2 cos{ujdt)) + (2x/C'e)^ -f 2dtx/(!/^) 

= ^ + ^(<^9 + #) +^('9t^')^0^- 

In the last line we used 2x1 Cy. = {dtr)/r and the 
first line of (30) to compute 2dtx/CY = {dtO)^ — 
Lol (1 + 'qcos{uJdt)) - (2x/Cs)2. 

The next transformation removes the cross term and 
rescales q and 0: 

[>4 = , (36) 

where F{t) = ^ log(r) satisfies dtF = x/Cs- Given 
[(0g -I- <70), 0] = —2ih(j) and [(0g -|- q4>),q\ = 2ihq, we 
compute 

lllqUi = e~'^^q = ^ 
r 

Ul(j)U4 = e^^4’ = r(j) 

-ihUldtUi = -7^(09 + #). (37) 


In this reference frame, we can interpret as the creation 
operator for the instantaneous eigenstates of HLc{t), 
whose energies are integer multiples of h{dt0){t). 


5. INTERACTION PICTURE 

Using the series of unitary transformations derived 
above, we now compute the operators 0 and q in the 
interaction picture with respect to Eq. (40). This serves 
two purposes. First, it will allow us calculate in the in¬ 
teraction picture the higher order terms i?(0) of equation 
(20), which we we dropped in order to the reach driven, 
quadratic Hamiltonian of equation (26). In the appendix 
we evaluate the size of these corrections in the rotated 
frame, and suggest a technique for minimizing their ef¬ 
fect. Second, going to the interaction picture will allow us 
to write the original interaction Hj = in terms 

of the effective annihilation operators of equation (39), 
allowing for a straightforward calculation of the coherent 
coupling strength in the next section. 

It is easier to first compute 0 in the interaction pic¬ 
ture; the first two transformations are the displacements 
involved in linearization, 

0 —>■ 0 -|- 0c , (41) 

while the third unitary Us = exp(zx02//i) has no effect 
on 0. The final unitary U 4 rescales 0 by the factor r, 
giving 

\/3w’^^“ + “t) + 0c(t), (42) 

where we have used (39) to express 0 in terms of the 
annihilation operator a. Finally, since in this rotated 
frame the circuit Hamiltonian is of the form H^c = 
h{dt9) (d’^d -I- 1/2), in the corresponding interaction pic¬ 
ture d —>■ de“*^. The final form of 0 in the interaction 
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picture is thus 


6. COUPLING STRENGTH 


</>c(t) + ^ 


(43) 

</>c(t) + ^ 


(44) 


where we have used / = re*^. Using this expression for 
the interaction picture flux operator, in the appendix we 
bound the effect of the higher order corrections R{(t)) to 
the linearized Hamiltonian of equation (24). A parameter 
of interest arising in this discussion is the relative size of 
the characteristic flux, which is set by 


In order to compute the effective coupling strength be¬ 
tween the circuit and ion motion, we must first analyze 
the Fourier spectrum of the characteristic function /. We 
begin by expressing the interaction Hamiltonian of (49) 
in terms of dimensionless factors, 

Hintiu) = ^^^^l{duf*a + dufa^) 

where we have used (()o = h/(2e) to get an expression in 
terms of the flux parameter 7 of equation (45). Here we 
have rewritten / in terms of the dimensionless param¬ 
eter u = uJdt/2, to better match the canonical form of 
Mathieu’s equation. 


Importantly, the final coupling strength is linearly pro¬ 
portional to this parameter, which we shall see limits the 
efficacy of our approach. 

The derivation of qmt is analogous to that of ipint- 
From the action of the four transformations Ui through 
C/ 4 , q takes the form 

qdt) +-q + C^{dtr)$ . (46) 

r 

Expressing these operators in terms of a of equation (39) 
then going to the rotating frame a —)■ gives 

q qc{t) + {g{t)a + g{t)*a)) , (47) 

where g{t) = (9tr — z^)e“*®. Using the Wronskian iden¬ 
tity lU = r'^dtO of equation (31), we see that g{t) = 
i9t(re“*®) = dtf*^ so the final form of q in the interaction 
picture is 

qiut{t) = qc{t) + {dtfa -k dtfa^) . (48) 


[dl + {A-2Q cos 2u)] / = 0, . (51) 

Using the substitution A = 4(a;o/wd)^, and Q = — 770 / 2 , 
this equation is equivalent to (27). By Floquet’s Theo¬ 
rem, / can be expressed as a quasi-periodic function 

/(zz) = e*^“^Cfce2*'=“, (52) 

k 

where the sum has period uq = tt corresponding to 
t = 2TTjuJd- Multiplying the cross terms of equation (50), 
we see that the coupling strength is proportional to the 
time-independent part of (the only term 

to survive under the RWA). This constant is exactly pro¬ 
portional to the Fourier component of / corresponding 
to the ion’s motional frequency, which is specified by the 
resonance condition 

{g, + 2k) = 2uji/ujd ■ (53) 

Accounting for the derivative duf* in this expression, the 
coupling strength is 

= ^^? 7 (M + 2A:)|cfc| = , (54) 


With this result, we may immediately compute the final 
ion-circuit interaction. 


Rintit) — C 
= ? 




dCs 

ezo 

C^d 


Qinti^) ^int(d) 


\[^{dtrd+djd^) 


(49) 




where Zi„t(t) = Zo{be -k with zq = ^ 

the characteristic displacement of the ion. Note that we 
have dropped the term qc{t)zint(t) by using the rotating 
wave approximation, since qc oscillates at characteristic 
frequency ujd^ uji. 


where in the second equality we used condition (53). 
Note that although H is not explicitly dependent on the 
driving amplitude 77 and frequency both 7 and are 
functions of these parameters since both are dependent 
on the characteristic function /. 

Before we compute the effective coupling strength be¬ 
tween the circuit and ion, we first rule out the presence of 
other resonances in equation (50). To do so we account 
for the ion’s micromotion, which can be derived from its 
equation of motion in a time-dependent potential. Using 
calculations analogous to those in Sections 4 and 5, one 

may show that Xmt(t) = Zq -kh.c.^, where 

h{t) is a periodic function whose fundamental frequency 
ujrf matches the trapping potential’s RF drive [10]. In 
the previous analysis we have approximated h{t) = 1 , 
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as this represents the largest Fourier coefficient of h{t), 
while other coefficients correspond to frequencies dis¬ 
placed from uJi by integer multiples of Typically 

ujrf ranges between 10 — 100 MHz, in contrast to the 
Wd « 1 GHz frequency spacing of the charge oscillations 
(equations (48) and (52)). Since the Fourier coefficients 
of both h{t) and f{t) are negligible at higher frequencies, 
the product Zint{t)qint(t) only has a resonance between a 
and W at frequency By the same reasoning, the two¬ 
mode squeezing terms {ab and a^b^) oscillate at frequency 
at least 2uji and may be dropped in the RWA. Likewise, 
no resonance exists between Zint{t) and the classical mo¬ 
tion, since it qdt) is a periodic function with fundamental 
frequency cod- Thus, the only remaining term after mak¬ 
ing the rotating wave approximation on is the 

desired interaction Hamiltonian, MlaW -\- h.c.. 

With equation (54) we can now evaluate the strength 
of the coupling for a driving strategy similar to that of 
Ref. [32]. Specifically, we set the drive frequency to be ap¬ 
proximately the difference between the circuit’s bare fre¬ 
quency and the ion’s motional frequency: ujd ^ too — 
Since the ion frequency oji <?; loq is much smaller than 
the LC frequency, the drive frequency tOd is nearly res¬ 
onant with the LC circuit. This means that even a rel¬ 
atively small drive amplitude leads to a mathematical 
instability in the system. Indeed, for rj sufficiently large, 
the characteristic exponent /r (equation (28)) describing 
the quasi-periodic function / gains an imaginary com¬ 
ponent, causing the interaction picture charge operator 
qint(t) to diverge over time. This instability is alleviated 
in the presence environmental dissipation, as the system 
dynamics are then damped, though for simplicity we ne¬ 
glect these effects in our analysis. As seen in Fig. 3, for 
Wo « Wd ^ the boundary between mathematically 
stable and unstable regions is set by ?7 < 2-\/wi/wo. For 
near-resonant driving to be mathematically stable, the 
parameter rj must therefore be perturbatively small. In 
this 77 <?; 1 limit, the characteristic exponent is equal 
to fi = 2wo/wd -I- 0{id) [39], so the solution of the fre¬ 
quency matching condition ^ + 2k = 2 wi/wd corresponds 
to A: = —I. Using the relations A = 4(wo/wd)^ « 4 and 
Q = —r/A/2 « —2rj, these parameters can be mapped 
to the canonical form of equation (51). The coefficient 
c_i = Q/4 + O(Q^) = — 77/2 -I- 0{rd) is known from 
standard perturbative expansions [39]. Substituting this 
value into equation (54), we conclude that in the per¬ 
turbative regime rj < 2 a/w^/wo, the coupling strength 
between LC and ion modes scales as 

^ (1 + 0 ( 77 )) 

~ ^ (’ + 0(v7M)) ■ (55) 

We note that this driving scheme may not be optimal. 
When LOd is set far from resonance and 77 ~ 1, the char¬ 
acteristic exponent /x can be stable and vary over a large 
range of values, allowing for the resonance condition (53) 



FIG. 3: Stability diagram of Mathieu’s equation in the reso¬ 
nant regime. The parameters Wi/wo and 77 map to the canon¬ 
ical variables of equation (51) as A = 4(1— Wi/wo)^ and 
Q = — 277(1 — Wi/wo)^. The shaded region corresponds to 
stable quasi-periodic solutions, i.e. those with characteris¬ 
tic exponent /x having no imaginary part. The red dot at 
uji = O.OOlwQ, 77 = 0.06 corresponds to the driving parameters 
used in Fig. 4. 


to be satisfied at stronger driving. Note that this also 
changes the Wronskian W (which changes 7 cx l/\/W) 
in a non-linear way, so a general analysis for the best 
driving parameters is difficult. Strong driving may be 
infeasible for experimental reasons, since it may require 
too large a Josephson energy (which is proportional to 
/3 > 77 ). 

Unfortunately, we find that the final effective coupling 
H is substantially weaker than in the proposal of Ref. [32]. 
In that work, the ion’s motional degree of freedom was 
also capacitively coupled to the circuit, but the char¬ 
acteristic equation of the circuit was instead modulated 
by a time varying capacitive element. The same driving 
scheme as above was used, leading to an effective cou¬ 
pling strength of the form, 

Klcap ~ ^ ^ 


Here Qo = \J jg the characteristic charge fluctua¬ 

tion in the driven circuit, and 77 the relative amplitude of 
the time-varying capacitance, C(t) = C'e(I -I- 77 sin(w£it)). 
For comparison, from equation (55) we obtain a coupling 
rate H ~ I Hz. Decoherence rates for these systems are 
expected to be on the order of kHz, while leading or¬ 
der corrections from linearization scale as fiw/zJywo/wi)^, 
both of which make this approach infeasible as currently 











described^. Note that we are making the same assump¬ 
tions as Ref. [32] about trap geometry = 0.25, d = 
25/im) and ion mass (m ~ 1.5 x 10“^®kg for ®Be+), 
as well as ion motional frequency {oji « 27r x IMHz, so 

zq = ~ nm). Finally, the parameter 7 ~ 1.3 is 

set by the circuit capacitance Cs = 46 fF and the Wron- 
skian IF « wq = 27r x 1 GHz. 

As seen from the above comparison, when the ion- 
circuit interaction is capacitive (proportional to charge), 
modulating the circuit inductance at frequency ujd ^ uji 
is significantly less effective than modulating the capac¬ 
itance. This results from the asymmetric dependence 
between charge and flux operators on the characteristic 
function /, as demonstrated by the interaction picture 
expressions for these operators (equations (44) and (48)). 
While the flux scales with /(t), charge is explicitly depen¬ 
dent on dtf{t). This means that the Fourier component 
of qint{t) oscillating at frequency Ui (and thus the time- 
independent component of Hintit)) picks up a factor of 
Wi/wd ^ 1 compared to (pintit)- The opposite is true for 
capacitance modulation (for which the roles of q and 0 
are reversed in transformations U3 and t/4), which ex¬ 
plains why it achieves a much larger effective coupling 
for a charge based interaction. 

7 . CONCLUSIONS 

We have studied a technique to coherently couple a su¬ 
perconducting circuit to the motional mode of a trapped 
ion by careful variation of the circuit’s inductance. We 


describe a means of tuning the inductance (an exter¬ 
nal magnetic flux with a Josephson Junction) and de¬ 
scribe an approximation mapping the circuit’s non-linear 
Hamiltonian to that of a driven harmonic oscillator. 
Notwithstanding corrections to this approximation, the 
mismatch between the inductive driving and capacitive 
interaction is the major reason why the resulting cou¬ 
pling strength is impractically small. We confirm this in 
a direct comparison with a capacitive modulation scheme 
for a specific driving strategy, though the general form of 
the coupling suggests our conclusions hold for a broader 
class of strategies as well. Indeed, equation (54) holds for 
any choice of drive parameters ujd,ri, and in fact can be 
applied more generally to any periodic modulation of in¬ 
ductance^. Conversely, our work suggests that for an in¬ 
ductive interaction (e.g. one based on mutual inductance 
between off-resonant circuits) an inductive modulation is 
the preferred approach. 
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8. APPENDIX - LINEARIZATION 
PROCEDURE 


The expression for <f)int{t) allows us to evaluate the 
higher order corrections to the linearized Hamiltonian of 
equation (26). Since these corrections arise after the first 
two transformations {{4>, q) ^ (pc, q + qc)), they can 
be expressed as Rint{t) = R{<Pint{t) — (pc)- Using relations 
(20) and Ej = PPl/L we obtain 


Rint (0 


k>3 


y^Pckit) 


(/EEl 

\y 2CsIU^o 




h 

LC^W 


k>3 


(a/* -baV) 


fc 1 

M ’ 


where we have defined the characteristic flux parameter 


1 I h 


(57) 
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1 ^ h 

Here Cfe(t) = EJ is defined according to 

equation (18), giving 

j3cz(t)= /3sin(^c(0/^o) = V cos^ 

(3c4{t)= I3cos{cl)c{t)/^o) =r]cos{uJdt) 

f3ck{t)= , 

where we have used (3 cos{(j)c{t)/(po) = r]cos{uJdt). 

To bound the error introduced by Rintit) we begin by 
considering only the k = 3 contribution. Using uJq = 
1 /LCy, this term can be written as 

- {'n/PYcos^{ujdt) (af* + aV)^ ■ (58) 

To analyze this term under the rotating wave approxi¬ 
mation, we note that our coupling scheme is premised 
on giving / a Fourier component matching the ion mo¬ 
tional frequency, Specifically, in the driving scheme 
described in the text, the parameters « wq — uji and 
r] < f3 1 are chosen such that / has the form 

f{t) = 1^1 + I - 6-“-^*) -h ...) , (59) 

where all other terms are of order 0(ff') and correspond 
to frequencies nujd (|n| > 2). Thus the only slowly rotat¬ 
ing term in equation (58) is of order hujQjr] (since IT « wq 
for these parameters) and rotates at frequency This 
slowly rotating part is the main contribution of (58) for 
evolutions over a time scale ~ l/w^, and we may compute 
its effective size over this timescale by using a second- 
order Magnus expansion [40, 41], 



FIG. 4: A plot of the periodic component of the characteristic 
function, p{u) = e“*^“/(M), in the weak driving regime. The 
blue, solid line and red, dashed lines represent the real and 
imaginary parts of p{u), respectively. Because the driving 
frequency LUd = 0.999 wo is nearly resonant with the bare LC 
frequency, even a small driving amplitude p = 0.06 causes 
significant deviations from the undriven case (for which pu = 
2'KUJot and p{u) = 1 ). 

the leading order contribution of Rint{t) by a factor of 
1/A^"^. For the resonance ratio uiojuji = 1000, an array of 
N ~ 100 junctions should sufhce to limit the effect of all 
higher order corrections. 


Rint{t) f>Wi ^7/3^^ ■ (60) 

Using a similar procedure, we can derive analogous esti¬ 
mates for the higher {k > 3) order terms as well. But 
because these terms pick up extra factors of 7 , equa¬ 
tion (60) represents the overall scaling of all higher order 
terms. This is true when the scale of |/| is on the order of 
^ 1 (as in Fig. 4), the characteristic flux is small (7 < 1), 
and there are only small fluctuations above the classical 
solution, (d^d) ~ 1 . 

One technique to reduce the effective size of correc¬ 
tions Rintit) is to replace the single Josephson junction 
with a linear array of these elements. In the simplest 
approximation[42], using a stack of N junctions we get 
that the Josephson Hamiltonian contribution is trans¬ 
formed as 


—Ej cos ( ^ I — 7 > —NEj cos ( I (61) 
\(t>oJ \NcPoJ 

In terms of the parameters in the definition of Rintit), 
this corresponds to the map /3 —>■ (3/N, ipo —>• NtpQ (or 
equivalently 7 —> j/N). From equation (60), this rescales 








